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For explaining certain features of atomic spectra, of atom-laser multiphoton 
transitions and of other atomic or molecular processes, a knowledge of matrix elements 
of polynomial radial functions between radial hydrogenic eigenstates is obviously of 
major importance (Moss 1972, Wong and Yah 1983a, De Lange and Raab 1991, 
Quiney et al 1997). Such importance comes about since the terms can be 
regarded as explicit expressions for interatomic potentials or as terms in a multipolar 
description of the interaction with the electromagnetic field. In that area a great deal 
of work has been done in nonrelativistic quantum mechanics (Nunez- Yepez et al 1997, 
1995 and the references therein). Nowadays, however, the sophisticated experimental 
techniques available make necessary the knowledge of such matrix elements between 
relativistic hydrogenic states. Though it is possible to compute directly such matrix 
elements (Kobus et al 1987, Wong and Yeh 1983b) the computations quickly become 
cumbersome thus techniques for evaluating any number of them starting from a 
few known ones are extremely convenient. For example, motivated by the need of 
evaluating the curvature-induced modifications of the hydrogen spectrum Bessis et al 
(1985) have contributed in this direction. In spite of this, however, there is not as 
yet any valid recurrence relation in relativistic quantum mechanics analogous to the 
general and useful ones customarily used in ordinary quantum mechanics. 

In this work we derive one such recurrence relation between relativistic radial 
hydrogenic states; this relation expresses the matrix elements of and Pr^ (where f3 
is the Dirac matrix) in terms of those of r^^^, (3r^^^, r^^"^, (3r^^^, r^^^, and (3r^~^ — as 
given in equations (17) and (20) below. The recurrence relations we obtain here can be 
regarded as a generalization of the non relativistic Blanchard rule (Nunez- Yepez et al 
1995, Blanchard 1974) to the bound hydrogenic radial eigenstates of Dirac relativistic 
quantum mechanics. Such relation can be found useful in, for example, different 
schemes of perturbation theory used in atomic calculations (Brack 1983, Dobrovolska 
and Tutik 1999) or for computing the effect of external fields on bound electrons 
(Wong and Yeh 1983a, b). Please be aware that we only intend to give the main 
results in this communication, for most of the details and other possible applications 
see Martinez-y-Romero et al (2000). 

The Dirac Hamiltonian, the Dirac equation for the stationary states of the 
hydrogen atom, and the corresponding energy eigenvalues, are 



Zaphc 



i7^^'(r) = E*(r), 



E = mc 



\ 



-1/2 



(1) 



(n-j- 1/2 + V(J + 1/2)2 J 



where a. and /? are standard 4x4 Dirac matrices in the Dirac representation (Bjorken 
and DreU 1964; Drake 1996, Ch 22 by Grant IP), Z is the atomic number, r 
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the relative distance between the electron and the nucleus, m the electron mass, 
c the speed of ligth, ap = e^/47reo^c the fine structure constant, e the electron 
charge, j = 1/2,3/2,5/2,... the total angular momentum quantum number, and 
n = 0, 1, 2, ... is the principal quantum number. The obvious rotational symmetry of 
Hd implies that the bound eigenstates of the hydrogen atom can be written as 



where = + 1/2, . . . , j — 1/2, j, is the 2;-projection of the total angular 

momentum quantum number, and XK,mA^^^) X-K,mA^j4>) are spinor spherical 
harmonics of opposite parity and k = —e{j + 1/2) is the eigenvalue of the operator 
A = P{1 + S-L) (Il = <T(8)/ = diag(<T, O"), a is the vector with Pauli matrices as 
components), which can be seen to commute with Hd (Drake 1996, Ch 22 by Grant 
I P). We found it convenient to define the quantum number e = (— l)-?+'~i/2 — ^.g. 
e equals +1 when / = j + 1/2 or equals —1 when I = j — 1/2 — and use it instead 
of parity for labeling the eigenstates, so I = j + e/2 and V = j — e/2, because the 
small component has the opposite parity to the big one (Martinez-y-Romero et al 
1998, 1999). Please notice that another notation for the spinor functions used in (2) 
y\,m, = XK,m, and y^[^^ = X-K,m, (Martinez-y-Romero et al 1998, Greiner 1991, 
Moss 1972). Writing the eigenfunctions in the form (2) completely solves the angular 
part of the problem, so we only need to cope with the radial part of it. Be aware also 
that, as we did in (2), we prefer to employ n, j and e rather than just n, and k to 
label the radial eigenfunctions. 

Let us begin establishing a relativistic version of the hypervirial result which 
is known to lead directly — though not straightforwardly because the computations 
are rather long — to the Blanchard recurrence relation in nonrelativistic quantum 
mechanics (Nuiiez-Yepez et al 1995). The radial Hamiltonian associated to (1) 
can be obtained using the squared total (orbital plus spin) angular momentum 
J2 = + <T-L + 3/4, and the fact that the operator L^, as applied on eigenstates 
of the form (2), is equivalent to the action of the operator j{j + 1) + e(3{j + 1/2) + 1/4 
upon the same states. Thus, using (a • r)(a • p) = (S • r)(Il • p) = r • p + iS • L, we 
get a ■ p = [Pr — ■^/?(ife + l/2)efe/r], so we finally obtain (Martinez-y-Romero et al 





2000) 



Hk = OLr Pr-ip—[jk + ^\ +/5m + F(r), 



(3) 



where we use units such that h = c = 1, V{r) is an arbitrary radial potential. A; is a 
label introduced for the sake of later convenience, and 
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ar 



a 





-1 



i d 

Pr = -T-r, 

r ar 



(4) 



(r) = 



^^rik,jk,^k (^) 



As it should be clear from the matrix expression for ar in the previous equation, we can 
use a 2x2 representation valid for the radial eigenstates given in (4) — where the now 
2x2 /3-matrix is just diag(l, —1) with numerical entries (Constantinescu and Magyari 
1971, p 382). Using this representation the purely radial Dirac equation reduces to 



m+{V{r)-E) Ze{j + l/2)/r-d/dr 
_Ze{j + l/2)/r + d/dr m- {V{r) - E) 



Fnjeir) 



= 0; 



(5) 



we have to point out, however, that using this representation is not strictly necessary 
and that all our results are representation independent. 

The key relationship needed for deriving the relativistic recurrence relation 
stems directly from equation (3) — compare with equation (2) and the nonrelativistic 
discussion that follows in (Nunez- Yepez et al 1995 pp. L526-L527). 

Let us first compute matrix elements of the radial function ^(r) = H2f{r) — f{r)Hi 
between radial eigenstates of the hydrogen atom in the energy basis, where /(r) is an 
arbitrary radial function and the Hi are the radial Hamiltonians appearing in (3). 
Evaluating such radial matrix elements we get 



{E2 - Ei){n2j2 e2|/(r)|niii ei) 

= {n2j2e2\H2f{r) - f{r)Hi\nijiei) 



where \nje) = 



F 



we have defined the quantities Af-^ = e2{2j2 + 1) ± ei(2ii + 1), and the expressions 
used for the {n2j2^2\fir)\nijiei) and the {n2j2^2\Pf{r)\nijiei) matrix elements are 



(n2j2e2|/(r)|nijiei)= / f{r)[F^{r)F;{r) + G,{r)G;{r)]dr, 

(7) 

(W2|/3/(r)|nijiei)= / f{r)[F^{r)F;{r)-G^{r)G;{r)]dr. 
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We then apply the result (6) to the function H2^{r) — ^{r)Hi, to obtain 



{E,-E,r{s,\f{r)\s,) = 



i^^l - ^mr) - fir) - ^nr)P - ^fir)^i+ (g) 
^f'ir)P + (^) m + Wm (^/'(r) + ^/?/(r)) ki); 

where the subscripts in the radial components, F and G, stand for the three quantum 
numbers n, j, e, and we have assumed 0. From now on, as we did in equation 

(8), we use the shorthand \si) = \niji €{) for the states of the system. It is to be noted 
that (8) can be regarded as the relativistic equivalent to the hypervirial obtained in 
(Nuiiez-Yepez et al 1995, equation (8), p. L526). 

However, in the Dirac case we are dealing with, at difference to what happens 
in the nonrelativistic case, equation (8) is not enough for deriving the recurrence 
relation between matrix elements of powers of r, we also need the following results 
(Martinez-y-Romero et al 2000): 

a) A second order iteration for certain non diagonal matrix elements, where we use 
for computing the elements, the radial function H2^{r) + ^{r)Hi 



{El-El){s2\f{r)\s,) = {s2\H2^ + ^H,\s,) = 

^21^21 p/\ o;„ / -P'/„^ 1 ^21 



(10) 



4^2 

b) the following matrix elements 

-i{E2 + E^){s2\arf{r)\s,) = -{s2\^ + /'(r) + 2fir)^ 

-^Pf{r) + 2iarV{r)f{r)\s^), 

c) the matrix elements of the radial function —i{H2arf{r) — arf{r)Hi), lead to 

-i{E2-Ei){s2\arf{r)\si) = 

A+ (11) 

{S2\ - fir) + -^Pfir) + 2iarf3mf{r)\s,). 

and, d) the matrix elements of the radial function H2Pf{r) + Pf{r)Hi, lead to 
(Marti'nez-y-Romero 2000) 
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(E, + E,){s2\mr)\si) ={s2\iM{r) - i^r^/lr) ^^^^ 

+ 2[m + /3F(r)]/(r)|si). 

Equations (6) and (8-12) are the basic equations for our problem. For more details 
on the computations please see Martinez-y-Romero et al (2000). 

Up to this point, our results are valid both for an arbitrary radial potential V{r) 
and for an arbitrary radial function /(r) but, to be specific, let us assume that V{r) 
is precisely the Coulomb potential, i.e. 



V{r) = ~, (13) 

where we have taken e^/47reo = 1. In this work we consider only the case of functions 
of the form /(r) = r^; so, putting /(r) = r'*'"^ in equation (6) and in equation 
(10), extracting the term with —iarA2i(3r^~^ from the former and the term with 
—2Xr^~^d/dr from the latter and substituting them into equation (9), but evaluated 
using /(r) = r^, we get 



{El - El){s,\r'\s,) = {s^l^^r'-' + ^(1 - X)Pr'-'+ 
Z [2iQ;^r^-^(l - A) - 2(^2 - Ei)r^-'^] - {E2 + Ei)Xiarr^-^\si) 



(14) 



We next use /(r) = r in equation (6), extract the term with 2iar/3mr ^ and 
substitute it into (11) but evaluated using /(r) = r^~^, to obtain the following result 



{E2 — El) 



4mA 



(S2|-(A-I)r 



^21 



A-2 



(s2|(-iQ;,r^~')|si) = 

^{E2-Eiy + ^/3r'-^\si); 
ZA21 ^ 



(15) 



then, we use /(r) = in equation (12) and follow the same procedure as above, to 
get instead 



{E2 — El) 



A 



{s2\{-iay-')\si) = (^2| - (A - iy-'+ 



A-2 



4m 



A 



21 /3^A-2 



AZm 
A 



(3r 



A-l 



2 m 

X 



(16) 



{E2 + Ei)(3r^\si). 



To continue, we evaluate equation (15) first using /(r) = ^ and then using 
r"^"^, next we extract the term — iQ!rA(£^2 + Ei)r^~^ from the former and the term 
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—2iZar{X — l)r^ ^ from the latter, to finally obtain, on substituting these extracted 
terms into (14), the recurrence relation 

3 3 

co(s2|r^|si) = ^Q(s2|r^->i) + ^rfi(s2|/3r^->i), (17) 

i=l i=2 

where the numbers Cj, z = 0, . . . 3 are given by 



_ {El-E!){E2-E,)A^, 
^ (^2 - El) A- - 4mA ' 

21 



Cl = 



2Z{E2 - EiYA^ 



{E2-E^)A^,-Am{X-ly ^^^^ 
e, = ^-A(A-l)- (^^ + ^^)^- 



C3 



4 ' ' {E2- Ei)A^i-4:mX' 

-2Z(A-1)(A-2)A^, 



{E2 - Ei)A^^ - 4m(A - 1) 
and the numbers di, i = 2 and 3, by 



4 = ^ 



(1_A)+ ^(^= + SOA£ 



{E2 - Ei)A2i - 4mA 
Z{X - l)A2lA+ 



(19) 



(E2-Ei)A2l-4m(A-i: 



Equation (17), together with the specific values for the and the da, can be 
regarded as the direct relativistic version of the Blanchard (1974) relation — compare 
with equations (10) and (11) in Nunez-Yepez et al (1995). These relations 
are valid inasmuch as the number + W2 + A + 1 is greater than zero, where 
Wi = >/(j7+^V2p^^^^^"F = V ~ Z'^cip- This condition amounts basically to the 
requirement that any integrand should be at least of the form 1/r^'^'^ with 7 > 
(Martinez-y-Romero et al 2000). 

The extra complication in the recurrence relation (17), that does not occur in the 
nonrelativistic case, is the explicit appearance of matrix elements of powers of r times 
the P matrix. In a way, this is just a matter of a sign change in an integral — as it 
should be obvious from equation (7). However, since it is not possible to avoid the f3 
dependence in (17), as it stands such equation does not really allow the computation 
of (s2|r^|si) in terms of the {s2\r^~^\si) i = 1,2, 3 as it undeniably happens in 
nonrelativistic quantum mechanics (Nunez-Yepez et al 1995), something additional 
is needed. 

For obtaining the extra information required, we only need to multiply equation 
(16) times (£^2 — -E'l) — 4mA/A2i and equation (15) times {E2 — Ei) —mA2i/ X, and then 
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substract the results, to finally obtain the lacking recurrence relation for the matrix 
elements involving the P matrix times r-powers, namely 



eo(s2|/9r>i) =bo{s2\r^\si) + 62(s2|r^"'|si) + ei{s2\(3r^-'\si) 

(20) 

+ e2{s2\/3r^-^\s^), 
where the numbers di and Cj i = 1, 2, 3 are given by 

bo =4A [(^2 - EiY - 4m'] , 
62 =(1-A) [(A,1)'-4A'], 

eo =2{E2 + Ei)[{E2 - Ei)A^^ - 4mA], ^21) 
ei =4Z[4mA - {E2 - Ei)A2i], 

62 =^[(A,-0'-4A']. 

The validity conditions of this recurrence relation is that the number {wi + W2 + X + l) 
be greater than zero, exactly as before. 

Equations (17) and (20) are together the useful recurrence relations for evaluating 
radial matrix elements between rclativistic radial hydrogenic states; in this sense they 
are thus the actual rclativistic generalization of Blanchard rule. 

The recurrence relation obtained above (17) pressupose Agi ^ 0, but, for studying 
the diagonal case, we must have ei = 62 and ji — j2 (in other words Ki — K2) which 
precisely imply that A21 = 0. In such instance we cannot apply that relation, we have 
to rederive the recurrence relation using as starting equations 

{E2-E,){nde\f{r)\n2je) = (mj e|(-ia,/'(r))|n2i e), (22) 



{E2-E,){nde\-iarf{r)\n2je) Hn,je\-f{r) + ^(3f{r)+ 



2r' ' (23) 
+ 2iarPmf{r)\n2je), 



and 



{E2 + E,){mje\/3f{r)\n2je) = (m j e| - ia,/3/'(r) 

+ 2{m + l3V{r))f{r)\n2je), 



(24) 



where, for avoiding possible misunderstandings, we have reverted to explicitly writing 
the quantum numbers in the states and A+ = 2e{2j + 1). To get the result we want, 
we evaluate equation (22) with /(r) = r^, then we put /(r) = r^~^ in equation (23), 
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and /(r) = in equation (24); thence, using essentiaUy a procedure similar to the 
used in the 7^ case outlined above, we finally obtain (Martinez-y-Romero et al 
2000) 

A+ 

[{E2 - E^f - 4m^] {mj e\r^\n2j e) ^ X^imj e\/3r^"^\n2j e) 

-4:m{nije\Pr^-^\n2je) - 2m{E2 + Ei){nij e\Pr^\n2j e) ^^^^ 

- X{X-l){nde\r^-^\n2je). 

This is the version of the recurrence relation (17) which is valid when A21 = 0. On 
the other hand, relation (20) can be directly written in the case A2;^ = 0, giving a 
equation entirely equivalent to (25); so, in the case where A21 = 0, we have relation 
(25) and relation (20) with the restriction A21 = 0. 

As happens with the non relativistic recurrence relation (Nunez- Yepez et al 1995), 
various of its particular cases have interest on their own; for example when A = 0, 
from (25) we directly obtain 

[(^2 - Elf - 4m2] 6n,n2 = -4^rn {riij e\-\n2j e) 

r (26) 

-2m {E2 + Ei) {nde\(3\n2je), 

where 5ij is a Kronecker delta. This relation (26) could be regarded as a relativistic 
version of the well-known Pasternak-Sternheimer (1962) rule of non relativistic 
quantum mechanics, which says that the expectation value between hydrogenic states 
of the 1/r^ potential, vanishes when the orbital angular momenta of the states 1 and 
2 coincide, i.e. when li = I2. We remark that in the relativistic case the expectation 
value of the 1/r potential (which corresponds to the square root of 1/r^) times 
does not necessarily vanish even when the total angular momenta of the two states 
coincide: ji = j2. This agrees with the known non relativistic fact that the Pasternack- 
Sternheimer rule is applicable to eigenfunctions of potentials whose energy eigenvalues 
depend only on the principal quantum number — which is not the case for the hydrogen 
atom in Dirac relativistic quantum mechanics. 

Furthermore, in the completely diagonal case {i.e. when ni = 77.2, ji — j2, and 
ei = 62; or just 77-1 = 712 and Ki — K2) we easily find ((O) stand for O's expectation 
value) 

m = - {(3V{r)) + E {I3)^z(^(^'^+E {13), (27) 
or, using (/?) = E/m (De Lange and Raab 1991), 

E^ = m {(3V{r)) +711^ = -mZ +m'^. (28) 



10 R P Martinez-y- Romero et al 

This last expression can be regarded as another version of the relativistic virial theorem 
in a form different from the one proposed by Kim (1967). Please notice that the first 
equalities in equations (27) and (28) above are valid for an arbitrary central potential 

V{r). 

Some other useful things that can be obtained from both the relativistic hypervirial 
[equation (8)] and the subsequent results [equations (9-12)], and from the relativistic 
recurrence relation presented in this contribution [equations (17) and (20)]; but such 
discussions as well as the many details left out from our presentation, will appear in a 
more detailed article (Martinez-y- Romero et al 2000). Let us pinpoint again that for 
our results to have meaning we need that the exponent in the recurrence relations, A, 
comply with (A + 1) > —W1—W2, where we assume that the Wi = a/ (ji + 1/2)^ — Z^aj^ 
are real numbers (Martinez-y-Romero et al 2000). This condition has a similar form 
to the one required by Blanchard (1974) in the non relativistic case. 

We have just learnt of the plans to test CPT and Lorentz invariance studying the 
ls-2s two photon transition in hydrogen, in which the recurrence relations proposed 
here can be of some help for the computations (Bluhm et al 2000). 

We acknowledge with thanks the collaboration of C Cisneros and the comments of J 
J Pena and V Gaftoi. Last but not least, we want to thank the friendly support of 
M Dochi, C F Quimo, U Kim, Q Chiornaya, C Sabi, M Sieriy, Ch Cori, S Mahui, R 
Sammi, M Becu, K Zura, and M Mati. This paper is dedicated to the memory of C 
Ch Ujaya, F Cucho, R Micifus, and B Ch Caro. 
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